Spiraling eutectic dendrites 



Tamas Pusztai,^Laszl6 Ratkai,^Attila Szallas,^ and Laszlo Granasy^'^ 

^Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary 
^BCAST, Brunei University, Uxbridge, Middlesex, UBS 3PH, United Kingdom 
(Dated: February 21, 2013) 

Eutectic dendrites forming in a model ternary system have been studied using the phase-field 
theory The eutectic and one-phase dendrites have similar forms, and the tip radius scales with 
the interface free energy as for one-phase dendrites. The steady- state eutectic patterns appearing 
on these two-phase dendrites include concentric rings, and single- to multiarm spirals, of which the 
fluctuations choose, a stochastic phenomenon characterized by a peaked probability distribution. 
The number of spiral arms correlates with tip radius and the kinetic anisotropy. 

PACS numbers: 68.70.+W, 81.10.Aj, 81.30.Fb 



Dynamic evolution of spiraling patterns observed in a 
range of physical, chemical, and biological systems in- 
cluding excitable media (such as cellular slime mold [T]), 
Belousov-Zhabotinsky reactions [2j, growth on atomi- 
cally flat interfaces [3 , binary eutectic systems ^ , and 
more recently, in ternary eutectic systems [5] has been 
exciting the fantasy of researchers for some time. While 
the details differ, diffusion and phase separation usually 
play a role. For example, aggregation of starving cells 
is controlled by propagating spiral waves of a chemo- 
attractant, often yielding multiarmed spiral patterns [1^. 
In binary eutectics, spiraling has been associated with 
specific anisotropy of the solid-solid interface [4], screw 
dislocations [6 , or osmotic flow driven fingering [7]. In 
turn, the newly discovered spiraling ternary eutectic den- 
drites emerge from the interplay of two-phase solidifica- 
tion with the Mullins-Sekerka-type diffusional instability 
caused by the third component [5]. This spiraling/helical 
structure has been identified as of interest for creating 
chiral metamaterials for optical applications via eutec- 
tic self-organization [8]. The complex microstructure of 
some ternary alloys is suspected to originate from eu- 
tectic dendrites [9]. Remarkably, multiarm spiraling has 
been reported experimentally in excitable media [1], in 
binary eutectics [T, and in Liesegang reactions [10], and 
theoretically in the FitzHugh-Nagumo model, in which 
the multiarm spirals form due to the attraction of single 
spirals [11]. It is yet unclear how general this behavior 
is, in particular whether multiarm spiraling is possible for 
ternary eutectic dendrites, and what governs the number 
of spiraling eutectic arms. 

In this paper, we show that a minimal phase-field 
model of ternary freezing is able to describe the spiraling 
ternary eutectic dendrites, and perform a detailed numer- 
ical study of this exotic growth mode. We demonstrate 
that the multiarm eutectic spiral patterns are robust, so 
they should be experimentally accessible, and that anal- 
ogously to the findings for Liesegang reactions [10 , the 
number of spirals results from an interplay of stochastic 
effects and the competition of nonlinear modes. 

The free energy of a minimal ternary generalization of 



the binary phase-field model (see e.g. Ref. [12]) reads as 
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where (j) G [0, 1], e^, ec, and w are constants, for the g{(j)) 
and p{(j)) functions see Ref. [12 , whereas c = (ci, 02,03), 
and the bulk liquid and solid phases are regarded as 
ternary ideal and regular solutions: 



/;,,(c) = ^c, [/i'^ + logci 



l,s 

ij ^i^j' 



(2) 



The equations of motion (EOMs) have been derived vari- 
ationally, yielding 
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for the concentration fields, where the q = 1 con- 
straint is automatically satisfied by our choice of the spe- 
cific values of 1 and —0.5 for the diagonal and off-diagonal 
elements of the 3x3 mobility matrix, M^. 

The dimensionless form of these EOMs has been solved 
parallel by finite differencing and explicit time stepping 
on a 3D grid, on a cluster of computers. Most of our 
simulations have been carried out in a directional solid- 
ification configuration. A temperature gradient was im- 
plemented by making the solid free energy temperature 
dependent as fs^i^z = fs^i — ^ {dfs,i/dz)^ where z is the 
coordinate along the direction of Vp sample pulling. Sam- 
ple pulling has been modeled by shifting the contents of 
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the arrays (j) and q by one voxel back in the z direction in 
each [{dx/dt)/vp]^^ time step, with boundary conditions 
^ = and c = Co on the high T and no flux boundary 
conditions on the low T side of the sample. To enable 
large enough simulation boxes in the direction of pulling 
(z), only Eqs. (4) have been solved far ahead of the so- 
lidification front, where (j) is sufficiently small (< 10~^). 
Since the anisotropy of the solid-liquid interface free en- 
ergy is weak for metals and the transparent system used 
in Ref. [5j, we have considered only kinetic anisotropy 
(of cubic symmetry, see [13^). Solidification has been 
started by including a slab of solid of length ^^(2/9), 
where Lz = Nzdx is the length in direction z |14] with a 
small hump at the center, whereas the initial composition 
of the solid has been (ci) = (C2) = 0.455 and C3 = 0.09 
realized by a random transversal {x — y plane) distri- 
bution of the two solid phases. We have opted for this 
starting condition, because simulations that follow the 
formation of the two-phase dendrites from fiuctuation- 
induced emergence of the Mullins-Sekerka instability of a 
fiat interface are prohibitively time consuming. We have, 
however, demonstrated the formation of surface undula- 
tions increasing with time, and found that the lower and 
upper unstable wavelengths are ~ 60 and ^ 320, with the 
fastest growing wavelength being around ~ 120. The par- 
allel computations of this study would have taken more 
than 600 years on a single CPU core. 

In our study, first we have explored the parame- 
ter space defined by composition, temperature gradi- 
ent, pulling velocity, interfacial free energy, and kinetic 
anisotropy, and optimized the conditions for growing two- 
phase steady-state dendritic structures. Typical condi- 
tions for such dendrites are summarized in Ref. [M] . For 
low pulling velocities {vp < 0.03) one finds a planar front 
with lamellar pattern, whereas at high enough pulling 
velocities {vp > 0.3), solidification takes place without 
apparent partitioning, though the dendritic structure is 
yet preserved. At even higher pulling velocities {vp > 0.8) 
partitionless growth with a fiat interface has been found. 
A typical two-phase dendrite [see Figs. 1(a) and 1(b)] has 
a rounded square-like transverse section in the x — y plane 
[Fig. 1(c)]; whereas in the fin directions (e.g., x — z)^ the 
longitudinal profile can be fitted with z = Zmax — \^\^ ^ 
where Zmax is the tip position, x is the distance from the 
axis of the dendrite, and u is ^ 1A9 ± 0.05 [Fig. 1(d)], 
somewhat lower than the u = 1.67 found experimentally 
for single-phase xenon dendrites [15]. (The perimeter of 
the dendrite has been determined by the contour line 
(j) = 0.5.) This steady- state shape has been achieved af- 
ter a transient composed of decaying oscillations of the 
tip-radius, tip-temperature, and the maximum of C3 at 
the tip [Fig. 1(e)]. To test further, how far the two- 
phase dendrites resemble the single-phase dendrites, we 
have varied the magnitude of the solid-liquid interface 
free energy (7s'l), via changing the free energy of the 
single component solid-liquid interface (7s'l,o) and eval- 
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FIG. 1: (color online) Two-phase spiraling dendrite grown 
under conditions given in Ref. ^J: (a) spiraling motif on the 
surface; (b) the helical structure formed by one of the solid 
phases, (c) contour lines showing the transverse sections at 
lOdx distances; (d) longitudinal section (dots), the best fit 
parabola (dashed), and the curve z = Zmax — I^T fitted to 
it (solid line); (e) maximum of C3 at the tip vs time; (f) tip 
radius vs. solid- liquid interface energy. 



uated the tip radius in the fin direction (Rup). The re- 
sults indicate Rup cx JsI^q^'^'^^ [see Fig. 1(f)], which is in 

~ 1/2 

a good agreement with Rup cx 7^^^ derived theoretically 
for single-phase dendrites [16 , and may indicate e.g. a 
negligible chemical contribution to jsl- Apparently, the 
shape of the two-phase dendrite is independent from the 
eutectic pattern forming the solid dendrite: target pat- 
terns, single- and multiple spiraling motifs do coexist on 
the same Rup vs. jsl,o curve [see Figs. 1(f) and 2]. 

The target pattern advances via alternating nucleation 
of the two solid phases, a mode expected to disappear at 
small undercoolings. It is more frequent for lower in- 
terface free energies, and becomes rare for 'ysL,o > 0.1. 
Besides the target pattern, a number of steady state "spi- 
raling" modes have been observed that display one to five 
arms [Fig. 2]. (On the surface of the dendrite spirals 
are realized by helical structures forming in the volume.) 
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FIG. 2: (color online) Eutectic patterns of two-phase den- 
drites, (a) Front view; (b) longitudinal and (c) transverse 
sections. From top to bottom, 7sl,o = 0.0295, 0.0147, 0.0354, 
0.0516, 0.0589, and 0.0810, respectively. The disorder in the 
tip region increases with increasing interfacial free energy. 



Owing to evident geometrical constraints, the steepness 
of the spirals increases with the number of the arms. The 
longitudinal sections are fairly similar for all modes, al- 
though weak systematic differences are observed. More 
characteristic are the front views and the transverse sec- 
tions: The individual modes (number of spirals) can 
clearly be distinguished [Figs. 2(a) and 2(c)]. We also 
find that once in the appropriate parameter domain, the 
spiraling two-phase dendrites are quite robust. We note 
that in the experiments, which were performed at low 
undercoolings, only the single-spiral mode has been ob- 
served so far [5|. The large number of spiral arms seen 
here probably follows from the large relative undercool- 
ing ^4j used in our simulations. 

The larger the tip radius, the larger is the number of 
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FIG. 3: (color online) Probability distribution of the steady 
state dendritic patterns from 20 different random initial two- 
phase patterns, (a) ^sl,o = 0.0354; (b) ^sl,o — 0.0589. 



spiraling arms [Figs. 1(f) and 2]. A closer inspection of 
the tip region reveals that no nucleation is needed for 
the single spiral mode, where the two-phase spirals orig- 
inate from a rotating "yin-yang" like motif at the tip of 
the eutectic dendrite. The modes with larger number 
of spirals become increasingly more complex, displaying 
alternating phase appearance at the tip. It is difficult 
to decide whether heterogeneous nucleation or growth 
around the phase occupying the tip is the mechanism by 
which the phases invade the tip. With the exception of 
the target pattern, where cones of one of the solid phases 
are not connected with other cones of the same phase, 
in spiraling modes the individual one-phase regions are 
interconnected with all volumes of that phase [Fig. 1(a)]. 

Although the number of spiraling arms {Narm) tends 
to increase with the solid-liquid interfacial free energy 
(7s'L,o) [Figs. 1(f) and 2], the steady-state pattern ap- 
pearing after the transient depends also on the initial ran- 
dom distribution of the two solid phases: Different steady 
state patterns are obtained starting from different (ran- 
dom) initial patterns. For example, at ^sl,q = 0.0354, 
these patterns include the target pattern, and single- to 
triple spirals (Fig. 3), showing a multiplicity of steady- 
state solutions for nominally the same conditions of which 
random initialization (representing here the cumulative 
effect of preceding compositional fluctuations) chooses. 
In other words, the thermal fluctuations decide which 
steady-state solutions are accessible for the system un- 
der a given set of operating parameters. Indeed, we have 
observed a similar stochastic behavior, when initiating 
growth with a chemically homogeneous solid slab, and 
adding noise representing fluctuations to the EOMs (a 
study inspired by Ref. [lO]). These features closely re- 
semble the helical Liesegang patterns, where the thermal 
fluctuations determine, which of the competing modes 
(helical, double helical, or non-helical) is realized [10]. 
These similarities raise the possibility of a universal be- 
havior for a class of multiarm spiral systems. Whether 
such a stochastic behavior prevails in other multiarm spi- 
ral systems requires further investigations. 

Next, we investigate, how the kinetic anisotropy influ- 
ences the number of spirals. We find that with decreasing 
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FIG. 4: (color online) Tip radius vs. kinetic anisotropy at 
jSL,o = 0.0295. The number of spiral arms tends to increase 
with decreasing anisotropy, although with some scattering. 



ber of spiral arms tends to increase with the tip radius 
or interface free energy, and decrease with the kinetic 
anisotropy. These findings are expected to instigate fur- 
ther experimental/theoretical studies on multiarm spiral 
systems and their stochastic nature. 
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FIG. 5: (color online) Spiraling dendrite formed at c\ — 0.355 
while gradually changing ci/c2 at C3 = 0.09 in the incoming 
liquid, (a) Front view; (b) longitudinal and (c) transverse 
sections. Note the majority phase channel at the center. 

anisotropy the tip radius increases followed by the num- 
ber of the spirals, which however shows some stochastic 
scattering (Fig. 4). We note, furthermore, that the ex- 
ponent V describing the shape of the dendrite tip changes 
from ~ 1.49 ±0.1 to ~ 2.1±0.1, varying between roughly 
the experimental value for xenon (1.67) and the rota- 
tional paraboloid (2.0) expected for isotropic case. This 
is combined with a change of the transverse section from 
a square of rounded corners to a circle. 

Finally, we explore how the two-phase pattern varies 
for off-eutectic compositions. Eutectic dendrites have 
been seen to form only close to the eutectic composi- 
tion. With a slow change of the liquid composition, one 
can move yet away from the eutectic composition, while 
retaining the spiraling structure (Fig. 5). Beyond a criti- 
cal deviation from the eutectic composition, the majority 
phase forms a channel at the centerline of the two-phase 
dendrite, a feature apparent in the experimental observa- 
tions 0. Remarkably, such patterns exist in the steady 
state, raising the possibility that their apparent lack of 
formation from random eutectic pattern, is only due to 
a long relaxation time, inaccessible for our simulations. 

Summarizing, we have shown that the ternary phase- 
field model naturally incorporates the spiral eutectic den- 
drites, and that such two-phase growth forms emerge 
between the domains of lamellar eutectic patterns and 
solute trapping. The two-phase dendrites behave anal- 
ogously to their single-phase counterparts, whereas the 
underlying eutectic pattern has little influence on the 
shape. A number of eutectic growth modes compete, 
including the target pattern, and single to multiple spi- 
rals, of which thermal fluctuations choose. The num- 
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